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ABSTRACT 

We present a 3 — 5 fim LBT/MMT adaptive optics imaging study of three Upper Scorpius stars 
with brown dwarf (BD) companions with very low-masses/mass ratios {Mbd < 25 Mj^p] Msn/Mi, 
« 1-2%), and wide separations (300-700 AU): GSC 06214, IRXS 1609, and HIP 78530. We combine 
these new thermal IR data with existing 1 — 4 fim and 24 fj,m photometry to constrain the properties 
of the BDs and identify evidence for circumprimary /secondary disks in these unusual systems. We 
confirm that GSC 06214B is surrounded by a disk, further showing this disk produces a broadband IR 
excess due to smaU dust near the dust subhmation radius. An unresolved 24 fim excess in the system 
may be explained by the contribution from this disk. IRXS 1609B exhibits no 3 — 4 fim excess, nor 
does its primary; however, the system as a whole has a modest 24 fim excess, which may come from 
warm dust around the primary and/or BD. Neither object in the HIP 78530 system exhibits near- 
to mid-IR excesses. We additionally find that the 1 — 4 fim colors of HIP 78530B match a spectral 
type of M3 ± 2, inconsistent with the M8 spectral type assigned based on its near-IR spectrum, 
indicating it may be a low-mass star rather than a BD. We present new upper limits on additional 
low-mass companions in the system (< 5 Mjup beyond 175 AU). Finally, we examine the utility of 
circumsecondary disks as probes of the formation histories of wide BD companions, finding that the 
presence of a disk may disfavor BD formation near the primary with subsequent outward scattering. 
Subject headings: brown dwarfs, circumstellar material, instrumentation: adaptive optics, open clus- 
ters and associations: individual (Upper Scorpius), stars: individual (GSC 06214- 
00210, IRXS 160929.1-210524, and HIP 78530) 



1. INTRODUCTION 

Nearly all stars are born with optically thick, gas- 
rich, dusty accretion disks that are thought to be analo- 
gous to the early solar nebula and comprise the building 
blocks of planet formation. Within 5 — 10 Myr, most 
of these disks disappear, having presumably converted 
their mass into larger, planetesimal-to-planet-mass bod- 
ies. The timescale for this process typically decreases 
with in creasing stellar mass ( Carpenter et al.|2006 Cur- 
et al. 2009). The youngest very low-mass stars and 
dwarfs (BD) are also typically surrounded by 



disks. Their disk lifetimes also tend to increase with de- 



creasing BD mass (Luhman et al 



J |Monin et aL||2010| ), thollg 
behavior may not be universa 



2008 Carpenter et al. 



1 some work suggests 



ehavior may not be universal ( Luhman et al. 2010| 



rie 
brown 



20091 

this 

(hereafter L2010). 

While many young, single BDs, like stars, are sur- 
rounded by optically-thick accretion disks, it is as yet 
unclear whether all BDs formed like stars, in particular 
those that are companions to stars in wide binary sys- 
tems. Between 10 — 30 % of BDs are believed to form 
in binary systems ( Burgasser et al. 2003 



2003 



Close et al. 



Joergens 2008 ) , a frequency significantly lower than 



Observations reported here were obtained at the LBT and 
MMT Observatories. MMT Observatory is a joint facility of 
the University of Arizona and the Smithsonian Institution. The 
LBT is an international collaboration among institutions in the 
United States, Italy and Germany. LBT Corporation partners 
are: The University of Arizona on behalf of the Arizona univer- 
sity system; Istituto Nazionale di Astrofisica, Italy; LBT Beteili- 
gungsgesellschaft, Germany, representing the Max-Planck Soci- 
ety, the Astrophysical Institute Potsdam, and Heidelberg Uni- 
versity; The Ohio State University, and The Research Corpora- 
tion, on behalf of The University of Notre Dame, University of 
Minnesota and University of Virginia. 



that for low-mass M stars and sola r-mass stars (> 30% 
and ^ 50%;|Fischer fc Marcy|ll992l [Janson et al.p012b 
Raghavan et al.||2010D. If wide BD companions formed 



2010). it wide iiU companions 
like other, more massive stellar binary companions, they 
should comprise the tail end of the binary mass func- 
tion (BMF). Previous studies of the stellar BMF have 
foun d that low mass ratios (M b /Ma < 10%) are rare 
(e.g. Reggiani fc Meyer] 2011). In contrast, recent ob- 
servations suggest the frequency of at least solar-type 
stars with very low-mass BD companions may be higher 
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than expected if such objects were drawn from a.n ex- 
trapolatio n of the stellar ]3MF ( Ireland et al.|[20lT Jan- 
son et aL]|2012a ). Additionally, some ot the lowest mass 
wide-separation companions have masses straddling or 
below the deuterium-burning limit often used to distin- 
guish between planets and brown dwarfs (e.g. Luhman 



e t al.|l2006[|Metchev fc Hillenbrand||2006| 

in nght of these complications, several other, more 
planet-like, formation mechanisms have been proposed 
for these wide, low-mass BD companions. In these sce- 
narios the BDs form from the primary's protoplanetary 
disk, either in situ or initially much closer to the pri- 
mary. The first possibility is that the BDs form as part 
of a system of massive objects orbiting the primary and 
are subsequently scattere d by one of these objects to 
their current location (e.g. Veras et al.||2009 ). However, 
a dearth of close BD companions (the so-called "brown 
dwarf desert" ) has b een shown observationally and sug - 
gested theoretically ( |Marcy fc Butler||2000[ |Bate[|200o| . 
Another possibility is planet-like formation m situ via 
direct col lapse withi n the primary star's protoplanetary 
disk (e.g. Boss|20lT |. However, the typical stellar proto- 
planetary oIsKraHIus is 20 — 200 AU, so BDs formed at 
separations > 200 At/ would require stellar hosts with 



unusually large disks (Andrews et al. 20091. Further 



more, the mass ratios and separations ot wide BD com- 
panions are distinct from those of RV-detected planets 



HR 8799, 13 Pic, Fomalhaut, and 


K And (Marois et al. 


2008al 120101 ILagrange et al.|2010l 


Kalaset al..20U8l|<Jur- 


rie et al.||2012 Carson et al.||2013 


), indicating that they 



mation mechanism (e.g. Kratter et al.|[20T0 ). 



High-contrast imaging observations ot these BD com- 
panions and constraints on their disk population may 
clarify their formation mechanism. If wide BD compan- 
ions at hundreds of AU separation form in situ as the tail 
end of the binary mass function (just like stellar compan- 
ions) , we expect that they have the same disk fraction as 
single BDs, since binarity at these wide separati ons ap- 
pears to not affect circumstellar disk evolutio n (Jensen 



et al. 
they 



1996[ [Pascucci et aL|[2008l [Kraus et al^|20l2[ ). If 



lave a more planet-like formation mechanism, they 



are still likely to be born with disks. In our own So 
lar System, this is evi denced by the systeru of coplanar 
moons a,round Jupiter ( Lunine fc Stevenson|19 82'; 'Canup' 
fc Ward 2002 ). However, if these BDs form initially much 
closer to the primary and are scattered to their current 
locations, their disk masses may not be great enough 
to survive to the present day, as we will demonstrate. 
Furthermore, an intermediate-mass perturber should be 
present close to the primary star and may be detectable 
in imaging data. 

The nearby Upper Scorpius OB Association (USco) 
provides a good testbed for these scenarios. Its youth 
(5-10 Myr, [Preibisch et al.||2002| [Pecaut et aT] 2012 



Song et al. 2012| an d large niimber ot members (iV ^, 



500, Preibisch et al 



2002) mean that a significant pop- 



ulation ot massive disks remains, allowing for a robust 



from early to late-type stars and BDs 


Carpenter et al. 


2006 


2009 


Luhman fc Mamajek 


2012 


. Its proximity 



(145 jpc) provid es good sensitivity at low disk and com 



2002[ ). About 25% of the M-type BDs in USco exhibit 
broadband mid-infrared ex cess emission from warm dust 
in circum(sub)stell ar disks (Scholz et al.|20 07l'Ria z et~al 
2009, 2012; Luhma FlT MamaJekpnT^r fDawson eFlT 
|2Pl3| ). Sensitivity limits have thus far precluded similar 
studies of less massive L-type members. 

Furthermore, USco includes three stars with re- 
cently discovered wide-separation, very low mass (ra- 
tio) BD companions: GSC 06214-00210 (GSC 06214), 
IRXS 160929.1-210524 (IRXS 1609), and HIP 78530; 



(Ireland et al. 



2011 Lafreniere et al. 2008 



2011j), one 

to r etain a 

hot circum(sub)stellar disk ( ^Bowler et al. [20TT ). The 



ot which, eye 06214B, has been discovered 



panion masses (de Zeeuw et al. 1999 Preibisch et al 



first two systems both are K-type primaries and L-type 
BDs, while the third is a B-type primary and M-type 
BD. These three systems constitute the lowest known 
mass ratio binary systems in USco: q « 1-2%, which is 
exceedingly rare in ste llar binaries, regardless o f the pri- 
mary star's mass (e.g. Reggiani fc Meyer||2011 ). In fact 
these are among the highest mass contrast non-planetary 
systems found to date, populating the region of parame- 
ter space between planetary systems and stellar binaries. 
By probing the circum(sub)stellar environments of these 
and other similar systems, it may be possible to discrim- 
inate between the various suggested formation histories; 
here we take the first steps towards that goal. 

To explore the circumstellar properties of these very 
low mass ratio pairs, we present and analyze new adap- 
tive optics-resolved near/thermal IR imaging in four 
bands [3.1 fim, 3.3 fim, L' and M'] as well as archival un- 
resolved Spitzer/MIPS-24 photometry. In each sys- 
tem we search for the presence of a hot circum(sub) stellar 
disk component around the BD, identifiable by broad- 
band infrared excess emission in our array of resolved 
1 — 5 fim data. This emission is due to small dust, with 
temperatures above several hundred Kelvin, located in 
the in nermost (hot) disk where accretion is likely to be 
active ( |Tuthill et al.|2001[ ). The MIPS- 24 /zm data probes 
the presence ot warm ('J e/ / ^ 100 K) dust in the system, 
although the large beam size of these observations pre- 
cludes us from determining whether warm dust emission 
originates from around the primary or the secondary (or 
both). 

Our study is organized as follows. Section [2] describes 
our sample in more detail, our observations, and our 
data reduction. In Section [3j we re-examine the com- 
panions' spectral types and assess their near and mid-IR 
excesses using these new data. With our high-resolution 
adaptive optics (AO) imaging we are also able to provide 
constraints on the presence of additional companions in 
the systems. We further examine the disk properties of 
GSC 06214B and the nature of HIP 78530B in Section |4) 
discuss the utility of disks as probes of scattering histo- 
ries in Section [51 and summarize our findings in Section 

m 

2. SAMPLE, OBSERVATIONS AND DATA REDUCTION 
2.1. Sample 

Here we summarize the basic system properties of the 
three newly-identified very low-mass (ratio) BD compan- 
ions that are the subject of our study: 

• GSC 06214A is a 0.9 M© K7 star. Its mem- 
bership is based on proper motion as well as on 
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the youth indicators: X-ray emission , Ha emis 



sion, and hthium absorption (Bowler et al. 2011 
and references therein). GS(J U621 4B orpits at 
2.20" (320 AU) projected separation. I reland et al 
( 2011[ ) confirmed its common proper motion over a 
3 yr baseline and assigned it a spectral type of LO 
and a m ass of 14±2 Mjur, bas e d on 5 Myr DUSTY 
models ( [Chabrier et al.|2000[ ). [Bowler et al.| ( |2011 1 
spectroscopicaliy confirmed the object, and further 
suggested it is still actively accreting from a proto- 
planetary disk based on a 4 nm IR excess and Pa/3 



emission. 



• IRXS 1609A is another 0.9 M© K7 star, with 
membership determined similarly to GSC 06214A 



(Lafreniere et al. 2008 and references therein) 
: KXS"TgO'gB~brbi^ts ar2.22" (3 2 AU ) projected 
separation. Lafreniere et al. (2010) measured 
its common proper motion over a 1 yr base- 
line and spectroscopicaliy identified it as a young 
L4 spectral type. Unlike for GSC 06214B, the 
near-IR (NIR) spectrum of this object shows no 



sign of active accretion. At 8 



+4 



M 



Jup 



(based 



on 5 Myr DUSTY models), this object is below 
the deuterium-burning limi t nomin ally separating 
brown dwarfs from planets ( [Spie gel et al. 2011), al- 
though it has a higher mass ratio and much larger 
separation than directly-i maged planets such as 



P Picb and HR 8799bcd e ([Lagrange et al.[[2010 
Marois et al.[[2008a| [2010[). Thus, tor simphcity we 



refer to this object as a BD 



• HIP 78530 A is a - 2.5 Mq B9 star ( Lafreniere 
et aL|2011 and references therein); its iJSco mem- 
bership is based on direct parallax measurements 
by th e Hipparcos survey (157± 13 pc, van Leeuwen 
2007). HIP 78530B was discovere d at a projected 
separa tion of 4.5", or ^ 700 AU, by Lafreniere et al. 



(2011). The authors confirmed its common proper 



motion over & 2 yr baseline and spectroscopicaliy 
classified the BD as a low-surface gravity M8. They 
used 5 Myr DUSTY models to estimate a mass of 
23 ± 3 Mjup. As with IRXS 1609B, this object's 
NIR spectrum shows no indication of active accre- 
tion. 

These three systems (Figure[T]) may probe two very dif- 
ferent epochs in stellar evolution for the primaries and 
different internal structures for the companions. The B 
star HIP 78530A is on or near the main sequence, while 
the K stars IRXS 1609A and GSC 06214A are still de- 
scending onto the main sequence. HIP 78530B is a BD 
whose mass is clearly above the deuterium-burning limit, 
while GSC 06214B is near, and IRXS 1609B is likely be- 
low, this limit (for an assumed age of 5 Myr) . 

2.2. Observations 

For resolved 3 — 5 /im observations of each system, we 
used the near/thermal- infrared imaging camera, Clio, at 
the 6.5 m MMT and/or the Large Binocular Telescope 
Interferometer (LBTI) at the 8.4 m LBT (in single-sided, 
non-interferometric mode). A summary of these obser- 
vations is presented in Table [T] At 24 /^m, we used 



o 



K7 



GSC 06214 



LO 14M 



Jup 



o 



IRXS 1609 



L4 



Jup 




HIP 7853C 



o 



M8 {M3) 23 (>75)M 



Jup 



200 400 600 
Projected separation [AU] 

Figure 1. Schematic diagram of the three USco system proper- 
ties (discussed further in Sections [3| and [4|| . Circle diameters are 
proportional to the log of the object's mass. The black line is the 
projected separation beyond which objects as luminous as the "B" 
components would have been detected in the literature or our data. 
Fill indicates disk properties: objects without 1 — 5 fim or 24 
excess are open circles; the dark fill of GSC 06214B signifies its 
1 — 5 excess and actively accreting disk; light shading refers to 
unresolved 24 /im excess in the system which may be around the 
primary and/or companion. 

archival Spitzer/MIPS data for all three systems. Be- 
cause the MIPS point-spread function (PSF) full-width 
half-maximum (FWHM) is 1.3 — 2.7 times the projected 
separations of each companion from its host star, none 
of the systems are resolved at 24 /im. 

2.2.1. LBT/LBTI/LMIRCam L' and M' Imaging of 
GSC 06214 and HIP 78530 

The L and M Infrared Camera (LMIRCam) is the 
3 — 5 iim c hannel of LBTI (Hinz et al 



2008 Skrutskie 



et al.[[2010). The field of view with an //1 5 beam is 11" 



square, with 10.7 m as pixels (calibrated Nov. 2011 by 
Rodigas et al. (2012)). We assume the instrument rota- 
tion angle is repeatable to 1°. Field distortion has not 
been measured observationally; ZEMAX models of the 
LMIRCam optics predict field distortions of < 2% over 
distances < 5". 

The LBT AO system employs a deformable secondary 
mirror in order to minimize the number of warm surfaces 
along the optical path. As a result, it benefits from in - 
creased sensitivity in the thermal IR (|Lloyd-Hart[|2000 ) . 
For system details see Esposito et al. (2011) (Ell). 'I'he 
LBTI AO hardware design is identical to that used in 
Ell, but the calibration is independent. We corrected 
200 aberration modes, rather than the 400 modes Ell 
used to achieve > 80% Strehl ratio (SR) in H band. 
At the time of these observations, the LMIRCam PSF 
was degraded due to residual instrumental astigmatism 
(which has since been mitigated). Nonetheless, because 
wavefront phase errors decrease with increasing wave- 
length, our observations approached the diffraction limit. 
During this run, we achieved a SR on LMIRCam of 
80% at M' on a CrB {R = 2.2) at airmass 1. 
Data recorded directly from the wavefront sensor itself, 
not subject to the LMIRCam instrumental astigmatism, 
indicate the AO system was delivering ~ 90% SR dur- 
ing these observations. The science observations were 
executed under more difficult conditions: low elevation 
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Table 1 

Summary of LMIRCam and Clio observations. 



Object Date Instrument Filter Wavelength Range (tim) Integration (min) 



GSC 06214 2011 May 26 LMIRCam A/' 4.60 - 4.97 30 

2011 Apr 22 Clio 3.1 fim 3.03- 3.16 45 

2011 Apr 22 Clio 3.3 fim 3.13 - 3.52 20 

IRXS 1609 2011 Apr 17 Clio 3.1 fim 3.03- 3.16 45 

2011 Apr 17 Clio 3.3 fim 3.13 - 3.52 30 

HIP 78530 2011 May 26 LMIRCam L' 3.42-4.12 15 



(high airniass), increased winds, and with fainter guide 
stars, and so were subject to poorer AO performance 
than the a CrB observations. Thus the sensitivities we 
derive in this paper are somewhat lower than those which 
can now be achieved with LBTI on more favorable tar- 



gets (e.g. 


Rodigas et al.|2012 


Skemer et al.|2012 1 


its recently improved AU ca 


ibrations. 



We observed GSC 06214 and HIP 78530 on UT 2011 
May 26. During this LBTI commissioning run, only the 
AO system on the right LBT aperture was fully oper- 
ational, so we observed in single aperture mode. The 
observing strategy for these, and all other LBT I obser 



vations, wa s Angular Differential Imaging ( ADI) ( Marois 



et al. 2006 ) : the instrument rotation angle was fixed, and 
the held was allowed to rotate on the detector over the 
course of the observations. We used a four position nod 
pattern for taking image data, with nods of 1.5 — 3". 
Due to the high sky background, we were limited to very 
short exposures: 87 ms at M' and 146 ms at L' . To limit 
data volume, 200 individual exposures were coadded per 
single saved image. We obtained 30 min of M' data on 
GSC 06214 and 15 min of V data on HIP 78530. The 
FWHM (average of major and minor axes of the astig- 
matic PSF) for the HIP 78530 (L') and GSC 06214 (Af) 
images were 0.12" and 0.18", respectively. 

2.2.2. MMT/Clio 3.1 fj,m and 3.3 fj,m Narrowband Imaging 
of GSC 06214 and IRXS 1609 



Cli o is also a narrow-field 3—5 ixm imager ( Sivanandam 
. At the time of these observations it was installed 



2006 



on MMT; it is now at the Magellan Clay telescope. The 
MMT AO system has an architecture similar to that of 
LBT AO, with a deforma ble secondary mirr or correct- 
ing 52 aberration modes (Brusa et al. 2004). In //15 
mode, the Clio plate scale is 29.9 ±0.1 mas pixel~^, and 
the field of view of the entire array is 15" x 30" (cali- 
brated July 2010 using the binary system HIP 88817). 
We assume the instrument rotation angle is repeatable 
to 0.1°. Plate scale distortion has not been measured, 
but ZEMAX models predict distortion of < 0.2% over 
5" separations. 

We observed IRXS 1609 and GSC 06214 on UT 2011 
April 17 and 22, respectively. In order to decrease time 
lost to detector readout, we used only a subarray with a 
field of view of 9" x 30". Integrations were 10 coadds of 
5 — 6 s at 3.1 iim and 40 coadds of 1 s at 3.3 /im. The 
data were taken in ADI mode in a 15" two-position nod, 
with 1 — 2" dithers between nods. These observations 
suffered from moderate wind from the South that in- 
duced vibrations in the telescope and caused smearing in 
some of the images. The resultant FWHM ranged from 
0.20" — 0.35", depending on the severity of the vibrations. 



On GSC 06214, we obtained approximately 45 min and 
20 min of on-source integration in the 3.1 fim and 3.3 fim 
narrowband filters, respectively. On IRXS 1609 we ob- 
tained 45 min of integration at 3.1 iim and 30 min at 
3.3 i-im. 

2.2.3. Spitzer/MIPS 24 fim Imaging 

Each system had been previously observed by MIPS. 
IRXS 1609 and HIP 78530 24 um photometry are pub- 



lished in Carpenter et al. ( 2009 ) , and we refer the reader 
to this publication tor a description of the observations. 
For GSC 06214, we analyzed archival MIPS data taken 
as part of the Spitzer Legacy program: Gould's Belt 
(PID 30574, PI: L. Allen), which aimed to provide a com- 
plete census of star-forming clouds in the solar neighbor- 
hood. GSC 06214 was imaged in two AORs: 20000768 
and 20001280, obtained with the fast scan rate (3 s per 
frame). The system was located near the end of the scan 
map; therefore, only 5 frames per pointing were obtained 
near the source position with an equivalent integration 
time of 15 s at 24 /irn. We detected the source only at 
24 iim due to the shallowness of the survey; we found no 
source in the corresponding 70 and 160 /Ltm maps. There- 
fore, we only discuss the 24 ^m result hereafter. Multiple 
frames were mosaicked based on the image WCS pointing 
information and subsampled by a factor of 2, resulting 
in a plate scale of 1.245" pixel~^. 

2.3. Image Processing 
2.3.1. LMIRCam and Clio Data 

Custom MATLAB scripts were used for reducing the 
LBTI and Clio data. Background, bias, and dark sub- 
traction were achieved by differencing images from ad- 
jacent nods. Bad pixels were fiagged using a nearest- 
neighbors median filter and replaced with that median. 
In LMIRCam, each of the channels had a small time- 
varying bias level offset as well as variable pattern noise. 
We subtracted the median value of each channel to re- 
move the offset. Then we subtracted a 2-D template of 
the channel read noise created from the median of all 
channels. For those calculations we masked out a ~ 2" 
square region containing the source to avoid biasing the 
median. Finally, we aligned the images according to their 
position angle, rotating so North was up. 

These observations did not achieve sufficient on-sky ro- 
tati on to synthe size and subtract the primary's PSF, as 
per 



Marois et al. 



^ ^ (2006) . However all secondaries were 

sufficiently separated trom their primaries (> 2") that 
contamination from the primary was negligible. To pro- 
duce the final image, we mean-combined our set of ro- 
tated images together, and 3 — 4cr outliers for each pixel 
in the image stack were clipped. 
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Aperture photometry was used to extract the differ- 
ential flux between primary and secondary. None of 
the objects were resolved, so we used aperture radii of 
~ 1 FWHM. The background was estimated using neigh- 
boring apertures of the same radius at the same separa- 
tion from the host star. No aperture correction was nec- 
essary, because the same aperture radius was used for 
both the primary and secondary. 

We next calculated sensitivity function of ra- 

dial separation from the primary star. Residual low- 
spatial frequency offsets were removed by subtracting 
a smoothed image (Gaussian FWHM=10A/Z3) from the 
science image. The resulting image was then smoothed 
with a disk the same diameter as the aperture used for 
photometry. The noise, cr, as a function of separation was 
calculated by measuring the standard deviation of 1 pixel 
wide annuli centered on the primary. From this we calcu- 
lated the signal to noise ratios (SNR) of our detections. 
The uncertainty in our photometry (in magnitudes) is 
given by 2.5 • (/n(10) • SNR)~-^. The noise counts as a 
function of radius were also converted to contrast (delta 
magnitudes between the primary star and the measured 
quantity, noise in this case) by comparing to the peak 
pixel value of the primary star in the smoothed image. 

The primary star photometry in our filters was esti- 
mated by interpolating between published photometry 
in adjacent filters (i.e.: between shorter and longer wave- 
length photometry), because no standard stars were ob- 
served. There are no strong spectral features in K7 or 
B9 stars between 2 fj.m and 5 fim which would bias this 
interpolation. The 2MASS and WISE surveys recorded 
unresolved photometry for each system. Additionally, 
Carpenter et al. (2006) conducted a Spitzer IRAC sur- 
vey which includeaTRXS 1609 and HIP 78530. As we 
will show, the BD contribution to the system flux at 
1 — 5 fj,m is neghgible, so we can treat 2MASS, WISE 
Wl and W2, and Spitzer IRAC photometry as measure- 
ments of the primary star only. Where available, we also 
used resolved photometry at additional wavelengths from 
previous work. The data used for interpolation are listed 
in Tables [2l[3j andg 

2.3.2. Spitzer MIPS Data 

We reduced the raw 24 yum data of GSC 06214 using 
the MIPS Data Analysis Tool (iGordon et al.||2005|) with 



post-pipeline processing to improve hat hel dmg and re- 
move instru mental artifacts as detailed in Engclbracht' 
et al. (20071. We extracted the photometry using PSF 
fitting with the IDL-based StarFinder program ( |Diolaiti| 



et al.||2000[ ). The input PSF was constructed using ob 
servecTstars and a smoothed STINYTIM model PSF; 
it has been tested to ensure the phot ometry results are 



consi stent with the MIPS calibration (Engelbracht et al. 
2007|. 



3. RESULTS AND ANALYSIS 

In this section, we present photometry and analysis 
for each of the three binary systems. We first measure 
photometry for both components, from both our newly 
obtained high-resolution near-IR imaging and archival 
unresolved Spitzer 24 fim imaging. Interstellar extinc- 
tion/reddening can mimic the red excess from a hot disk 
component; the average extinction towards USco is small 



[Ay < 1), but not negligible, and we estimate it as nec- 
essary. We re-examine the companions' spectral types 
and constrain their near-IR (hot dust) excesses based on 
their placement in color-color diagrams compared to sim- 
ilarly young low-mass objects. We use unresolved Spitzer 
24 fj,m data to probe the warm dust in the system as a 
whole (both A and B components together). Finally, we 
quantify our sensitivity to additional objects in the sys- 
tem at intermediate separations from the primary. 

3.1. GSC 06214 
3.1.1. Photometry and Astrometry 

The final Clio images at 3.1 /im and 3.3 /xm are shown 
in Figure [2j At 3.1 /xm we measured a contrast be- 
tween "A" and "B" of 5.26 magnitudes, and at 3.3 jjLm 
we measured a contrast of 5.30 magni tudes. Using the 
procedure described in Section |2.3.1 wc inferred the 
3.1 /im and 3.3 /im magnitudes of the primary to be 
9.11 ± 0.05 and 9.10 ± 0.05, respectively. The 5cr sensi- 
tivity limits for these two filters were magnitudes 15.82 
and 16.14, yielding signal to noise ratios of 18 and 23 for 
the detections in each wavelength. From this we derive, 
^.1]b = 14.41 ± 0.08 and [3.3]b = 14.40 ± 0.07 (Table 

We calculated a separation of 2.17±0.02"and a position 
angle of 175.9±0.1°, co nsistent at the 1.5a level with 
astrometry measured by Ireland et al. (2011). Because 
the companion was not detected in individual frames, 
we estimated the uncertainty by comparing the values 
calculated from the final images at each wavelength. Our 
quoted uncertainty accounts for this comparison as well 
as our assumed uncertainties in calibration and distortion 



(see Sectio n |2. 2. 2[ ). 

The final LMIRCam Al' image is also presented in Fig- 
ure[2j The BD was detected at a contrast of 4.65 magni- 
tudes. For the primary, we estimated Af^ = 9.10±0.05. 
Based on the ba sensitivity limit of 13.39 mag, the 
BD detection has SNR = 3.6, and therefore M'^ = 
13.75±0.3. Even though this SNR is lower than that 
nominally adopted to identify can didate comp anions in 
high-contrast imaging data (5 — 6cr; Marois et al. 2008a), 
we also detect point sources at the same location in the 
3.1 /im and 3.3 /im images, justifying a lower detec- 
tion threshold. The bright spot ^ 1" to the east of the 
primary does not appear in the more sensitive 3.1 and 
3.3 /im images and is almost certainly an artifact. 

The 24 flux density of the system is 3.3 ± 0.3 mJy 
(8.34 ± 0.1 mag)] the error is dominated by shot noise. 

3.1.2. Extinction Estimate 

GSC 06214A did not have a published extinction value, 
so we performed our own estimate. We found Ay < 0.6 
using the following prescription. We compared the pri 



et all 2011 
' — . rrt' — 1 1 II ' 


) with 


Pickles 


(1998 



y their 800 — 900 nm a verages. Interste llar extinc 



tion curves with Ry = 3.1 (Cardelli et al. [1989 ) were 
used. Based on fltting, both K7V with Ay = and 
K5V with Ay = 0.6 provide acceptable fits to the data. 
K4V and Ay = 0.85 is a worse fit to both the continuum 
and spectral features. In order to make the most conser- 
vative estimate of the IR excess in the system, we adopt 
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Figure 2. GSC 06214 Clio and LBTI images. Left: Clio 3.1 fim. Center: Clio 3.3 ^tm. Right: LBTI A/'. The companion, circled, is 2.2" 
to the south of the primary. For display purposes, all images have been smoothed with a Gaussian kernel with FWHM=1A/D ~ 0.1". In 
the M' image, the bright spot to the left of the primary is an artifact; the dark spots to the left and bottom right are residuals from nod 
subtraction. 



Table 2 

GSC 06214-00210 system properties and 
apparent photometry. 



Property 


GSC 06214A 


GSC 06214B 


Distance [pc] ^ 


145±20 


Spectral type ^ 


K7±0.5 


LOitl 


Mass ^ 


0.9±0.1 Mq 


14±2 Mj^p 


T.J J [K] ^ 


4200±150 


2200±100 


logiL/LQ) b 


-0.42±0.08 


-3.1±0.1 


Separation ["] 
PA [°] 


2.17±0.02 


175.9±0.1 


J 


9.998±0.027 


16.24±0.04 " 


H 


9.342±0.024 


15.55±0.04 «^ 


K 


9.152±0.021 


14.95±0.05 <= 


3.1 fim 


9.11±0.05 


14.41±0.08 


3.3 fim 


9.10±0.05 


14.40±0.07 


3.4 ^lm f 


9.083±0.022 


L' <= 


9.10±0.05 


13.75±0.07 


4.6 fim ' 


9.107±0.020 


M' " 


9.10±0.05 


13.75±0.3 


12 fim f 


8.964±0.032 


24 fim 


8.34±0.1 



^ Mean cluster distance | |de Zeeuw et aL]|1999| 
[Preibisch et al..2002l . 



Bowler et al. 
i'nis work. 



120111 



2MASS J/H/Kc! su rvey jSkrutskie et al.|2006 i. 
'' Ilreland et al] | |2011| . MKU JJH/K'/L'. 
f System. WISE survey jWright et al.j2010 l. 

the maximum allowed extinction value: Ay = 0.6 mag. 

3.1.3. Spectral Classification and Evidence for a Near-IR 

Excess 

We placed GSC 06214 in infrared color-color dia- 
grams, compared to several observationa l brown dwarf 
and p lanet samples: field brown dwarfs 
20101, very young substellar objec ts (L20 



Leggett et al. 
0), the young 
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Figure 3. GSC 06214A flux vs. wavelength for best fit spectral 
type and extinction. Fx is offse t vertically for clarity. Gray is 
the stellar spectrum (Bowler et al. 2011), red is a K7V with no 
extinction, blue is a K5V with Ay = 0.6 mag, and green is a K4V 
with Ay = 0.85 mag. We find Ay < 0.6 mag for the primary. 

Chameleon, and TW Hya, so it provides the most appro- 
priate comparison sample for our similarly young brown 
dwarfs. 

In order to conduct a homogeneous analysis, we trans- 
formed all previously published J, H, and K BD pho- 
tometry to the MKO J/H/K system. For 2MASS 
J/H/Ks — )■ MKO J/H/K , we u sed the transformations 
pubhshed in [Leggett et al.| (|2006[). For MK O K' ^ K of 
stars, we uswTWainscoat & Cowie (1992); however, we 
could not apply this transformation to the BD, because 
it is only calibrated for H — K < 0.4. Therefore, we cal- 
culated the K' ^ K transformations by integrating over 
DUSTY evolutionary model sp ectral energy distributions 



(SEDs) (Chabrier et al. 2000). Wc took the SED that 



planetary-mass object 2M 1207b ([Chauv in et al. 2004; most closely matched the color of the BD in our filters of 
Mohanty et al. 2007; Skcme r et al.| |20lT; T3arman et al. interest (in this case, H — K). We used the same proce- 

' ' ' dure for all color transformations for which empirical re- 

lations were not available. We calculated K = 14.75 mag 
for GSC 06214B. We then dereddened aU GSC 06214B 
photometry, using the extinction value for the primary 




was derived from dereddened observations of Taurus 
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{Ay = 0.6) discussed in Section 3.1.2 We presumed the 
extinction toward the primary and secondary was the 
same, a nominally valid assumpt ion, though counterex - 
amples do exist (e.g. GG Tau, [White &: Ghezl|2001[) . 
Figure H shows the color-color diagrams H — L' / J — H 
and H - M'/J - H. 

GSC 06214B was previously clas sified using JJH/K 
spectra as a spectral type of L0±1 ( Bowler et al.||20lT ). 
Its J — H color is consistent with an L2U1U spectral type 
of M9. The excesses at L', and M' compared to an 
L2010 M9 dwarf photosphere are 0.15±0.05, 0.35±0.07, 
and 0.4 ± 0.3 mag, respectively, accounting for photo- 
metric errors only. Fortuitously, the measured excesses 
are similar even if we use Ay = 0, because reddening 
moves the object up and right in color-color space par- 
allel to the line connecting L2010 M9 and LO spectral 
types. Therefore, our conclusions about the color excess 
for this object are independent of the extinction estimate, 
over the allowed range of extinction values. 



3.1.4. 24 /im Flux 

The GSC 06214 system was imaged with Spitzer/MIPS 
at 24 /XTO. The PSF is broader than the host-companion 
separation, so it is not possible to measure the primary 
and secondary objects' emission independently. How- 
ever, the BD's photosphere does not contribute signifi- 
cantly at this wavelength. Whether based on a Rayleigh- 
Jeans extrapolation from short wavelength data or using 
DUSTY model SEDs, the BD's photospheric emission is 
expected to be less than 1% of the observed 24 iim flux. 

We compared the MIPS photometry to the 24 jim flux 
predicted for the "A" component by an empirical poly- 
nomial relation between V — Kg and 24 for main 
sequence stars ( Urban et al.||2012 ). The 24 jjLm flux pre- 
diction is robust to order 0.01 mag between early B and 
early M spectral types. GSC 06214A is likely in the final 
stages of contraction onto the main sequence. However, 
the variation 'm.V — Ks color, and therefore in inferred 
24 jim flux, between this object and a main sequence 
star, is less than the combined observational uncertain- 
ties at V, Ks, and 24 ^m. 

From this comparison, we find evidence of a 24 fj,m 
excess in the system. The observed system flux is 
3.3 ± 0.3 mJy, while the predicted flux from "A" is 
1.73±0.15 mJy. This translates to an excess of 91±24% 
above the expected photospheric emission from "A." The 
uncertainty is dominated by the low signal to noise of the 
24 ^m data. Our detection corroborates the 2a excess at 
WISE 22 fim found by 



Table 3 

IRXS 160929.1-210524 system properties and apparent 



(2011). Therefore, 



Bowler et al. 

one or both members oi the system must retain a warm 
excess. 

3.1.5. Constraints on additional objects 

No additional objects are de tected in our observations. 
As discussed in Section 3.1.1 the candidate at 1" sep- 
aration in M' is likely an artifact. In our collection of 
images, we could have detected additional companions 
of mass > 20 Mjup at separations > 150 AU from the 



photometry. 


Property 


IRXS 1609A IRXS 1609B 




145it20 


Spectral type ^ 




Mass ^ 


0.85l[];?[j Mq SttMj^p 


rp b 


40601^0° 1800l?og 


log{L/LQ) b 


-0.37±0.15 -3.55±0.2 


Separation ["] 
PA [°] 


2.15±0.03 


28.3±0.4 


J 


9.820±0.027 17.90±0.12 


H 


9.121±0.023 16.87±0.07 


K 


8.916±0.021 d 16.19±0.05 


3.1 fim 


8.80±0.05 15.65±0.21 


3.3 fim 


8.78±0.05 15.2±0.16 


3.4 fim s 


8.767±0.023 


L' <=>f 


8.73±0.05 14.8±0.3 


4.5 fim ^ 


8.80±0.01 


4.6 fim s 


8.779±0.021 


8.0 fim 


8.735±0.005 


12 fim s 


8.715±0.027 


24 fim ' 


8.43±0.01 



^ Mean cluster distance l |de Zeeuw et al.|1999| |Preibischl 
let al.|2002|l. 

|Lafreniere et al.| pOOs] !. 

1 his work. 

2MASS J/H/Ks survey i Skrutskie et al.|2006 l. 
" MKO J/H/K'/L '. 



Lafreniere et al. 



,1 2010 1 

System. WISE survey jWright et al.|2010t. 
^ System. Spitzer IRAC jCarpenter et al.|2006 1. 
' System. Spitzer MIPS ( [Carpenter et al.|2009l l. 

panion to "B" (i.e. "Bb") if it was > 20 AU projected 
separation from "B." 

3.2. IRXS 1609 
3.2.1. Photometry and Astrometry 

Figure [5] shows the final IRXS 1609 narrowband im- 
ages. We estimated the primary's 3.1 /zm and 3.3 ^.m 
magnitudes to be 8.80 ± 0.05 and 8.78 ± 0.05, respec- 
tively. At 3.1 /iTO we found a contrast between "A" and 
"B" of 6.85 mag. The 5(t background sensitivity limit 
of 15.71 mag yielded SNR = 5.3 for the detection, thus 
[3.1]b = 15.65 ± 0.21. At 3.3 /im we found 6.4 mag con- 
trast, a sensitivity limit of 16.14 mag, SNR = 7.5, and 
[3.3]b = 15.2 ± 0.16. (See Table |3]). 

We measured a separation of 2! 15" ± 0.03" at PA — 
28.3° ± 0.4°, 
the 2(7 level. 



consistent with Lafreniere et al. (2010) at 
These values represent the average of the 
astrometry calculated at 3.1 //to and 3.3 /im. Errors 
account for measurement uncertainties as well as our as- 
sum ed cah bration and distortion uncertainties (see Sec- 
tion [2^21). 



3.2.2. Spectral Classification and Near-IR Colors 
There is no exti nction towards the IRXS 1609 system 



(Carpenter et al. 2009, ). We transformed the K' com- 



primary, with a limit of ~ 10— 12 Mj^p beyond 200 AU. panion photometry to yield K = 15.99 ± 0.1, following 

the procedure described in Section [3.1.3[ IRXS 1609B is 
shown in the color-color diagrams in Figure |4] 

Unfortunately, the L2010 sample does not extend past 
spectral type LO; we therefore compared IRXS 1609B 



Previous aperture masking and direct imaging work have 
placed deeper limi ts of 35/15/5 Mj^„ at 10/150/300 AU 
from the primary ( Kraus et al.|2008 Ireland e t al.|20lT ) 
We would have resolved an equally bright bmary com 
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Figure 4. H — L' vs. J — H (left) and H — M' vs. J — H (right) color-color diagrams of dere ddened USco compan ions and comparison 
samples. Open green squares, brown diamonds, and gray circles: field M, L, and T dwarfs from'Leggett et al.' fZOlO"); filled green squares 
and brown diamond: voung MO — LP dwarfs from L2010; filled blue circles: young planets HR Svyybcd (Uurric et al. 20 1l[ [Marois et al.| 
|2008a[ [20T0I [Skemer et al.|2012| |Galiclier et al.|2011i l; filled green circle: young planet 2M 1207b ( jMohanty et al.,20U7| >. 





Figure 5. IRXS 1609 Clio images. Left: 3.1 fim. Right: 3.3 fim. For display purposes, both are smoothed with a Gaussian kernel with 
FWHM=1A/D ~ 0.1". The companion, circled, is 2.2" to the NE. 



broadband colors to DUSTY model SEDs with tem- 
peratures of 1700if and ISOOX. A surface gravity of 
^og{g) = 4.0 is predicted by the DUSTY models for these 
effective temperatures at 5 Myr. We scaled the SEDs so 
their integrated J band flux matched the observed flux. 
Both SEDs provide an acceptable fit to the dat a (Figure 



(120101) effective 



[g]), consistent with the Lafreniere et al. ^ 

temperature estimate. Our new photometry confirms the 
lack of near-IR excess. 

3.2.3. 24 fim Flux 

Following the procedure used in Section [3.1.4[ we find 
evidence for a moderate 24 fim excess in the system. The 
observed flux is 3.06 ± 0.04 mJy, and the predicted flux 
from "A" is 2.15 ±0.07 mJy, an excess of 42 ±6% above 
the expected stellar photospheric emission. The error is 
dominated by the uncertainty in the literature V mag- 
nitude. As with GSC 06214, the predicted contribution 
to the excess from the BD's photosphere is < 1%. Nei- 
ther "A" nor "B" exhibits a near-infrared excess, and 



the components have not been resolved at wavelengths 
longer than L' . Thus, though we can detect the presence 
of a warm excess in the system, we cannot determine 
whether it is around the primary, BD, or both. 

3.2.4. Constraints on additional objects 

We do not detect any additional objects in the 
IRXS 1609 system. Our background limit at 3.3 
was 15.6 mag, or ~ 5 Mjup, extending out to 500 AU. 
We do not improve upon previous detectio n limits of 
8/1 Mj„j, a t 50/450 AU stellocentric distance (Lafreniere 
et al. 2010 ). We could have detected an equally lummous 
binary companion to "B" (i.e. "Bb") at a separation 
> 25 AU from "B." 

3.3. HIP 78530 
3.3.1. Photometry and Astrometry 

Figure [t] shows the final L' image of HIP 78530, with 
the BD to the southeast. For the primary, we inferred 
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Figure 6. IRXS 1609B photometry and DUSTY SED comparison 
witli Ay = 0. Tlie green squares are ttie literature pliotometry and 
the blue diamonds are this work. Spectra are DUSTY 1800K and 
1700K model photospheres with log(g) = 4.0. The filled circles 
represent the integrated flux of each spectrum in that bandpass. 



Table 4 

HIP 78530 system properties and apparent 
photometry. 



Property 



HIP 78530A HIP 78530B 



Distance [pc] 

Spectral type ^ 

Mass ^ 

Teff [K] 1^ 

log{L/LQ) b 

Separation ["] 

PA [°] 

J 

H 

Ks 

3.4 ^lm " 
L' '^■^ 

4.5 fim ^ 

4.6 iim " 
8.0 iim s 
12 nm 
24 urn ^ 



157±13 
B9V M8±l 
~2.5 Mq 23±3 Mj„p 

~10500 2800±200 
-2.55±0.13 
4.54±0.09 
140.7±1 
6.928±0.021 d 15.06±0.05 
6.946±0.029 14.39±0.04 
6.903±0.020 14.17±0.04 

6.842±0.034 
6.87±0.05 13.80±0.06 
6.88±0.01 
6.890±0.019 
6.948±0.008 
6.924±0.017 
6.846±0.01 



Hipparcos catalog l|van Leeuwen|2007^ 



Lafreniere et al. 
i nis work. 



120111 



2MASS J/H/Ks survey fSkrutskic et a Lj2006 1 

lOf. 



System. WISE survey ( , Wrigh t ct al. 201 
f MKO L'. 



^ System. Spitzer IRAC I Cargenter^et^^ 
^ System. Spitzer MIPS i Carpenter et al. 



20061 



20091 



L'^ = 6.87±0.05. HIP 78530B was detected at a contrast 
between "A" and "B" of 6.93 magnitudes (Tabled. The 
5a sensitivity limit was L' = 16.21, giving SNu — 47 
and L'g = 13.80±0.06. Our measured separation and PA 
are 4.54" ± 0.09" and 140.7° ± 1°, includ ing pl ate scale 
and distortion uncertainties (see Section 2.2.11, consis- 
tent with previous observations. 

3.3.2. Spectral Classification and Near-IR Colors 
The extinction towards HIP 78530 is Ay = 0.5 (iCar- 



penter et al. [2009), and we dereddened the photometry 
using this value. We calculated a K magnitude for the 



Figure 7. HIP 78530 LMIRCam U image, smoothed for display 
purposes with a Gaussian kernel with FWHM=1A/D ~ 0.1". The 
companion, circled, is 4.5" to the SE. 

BD of 14. 12 ± .06 following the prescription described 
in Section [3.1.31 

We find the spectral types derived for HIP 78530B from 
broadband 1 — 4 /xm colors and from near-IR spectra are 
not consistent. HIP 78530B was classified u sing J/H/K 
spectra as M8 ± 1 (iLafreniere et al. 2011), though the 
authors noted the object was ~ 0.2 mag too blue at K 
for this spectral type. We find that its H — L' color is 
also ~ 0.5 mag too blue compared to an L2010 M8-type 
BD (Figure H left panel) . There is a similar mismatch 
in both H —~K and K — L' colors (Figure [sj left panel). 
With the addition of our new L' data, we find the object's 
J through L' colors are not consistent with an M8, but 
instead with an L2010 M3±2. Intriguingly, the object's 
absolute magnitude is consistent wi th other USco M8- 

fpe BD ( Luhman fc Mamajek|20l"2 ) as shown in Figure 
(right panel). We discuss the implications of this dis- 
crepancy further in Section |4.2| Regardless of whether 
we adopt spectral class M8 or M3, we find no evidence 
for K or L' excesses. 

3.3.3. 24 /im Flux 

The observed and predicted fiuxes for the HIP 78530 
system are 13.10 ± 0.12 mJy and 12.0 ± 0.5 mJy. The 
excess is 9 ± 4.5%. At 2a, we do not consider it statisti- 
cally significant. This, combined with a lack of 1 — 5 iim 
excess around both "A" and "B," suggests that neither 
object retains a massive warm disk. 

3.3.4. Constraints on additional objects 

No additional objects are detected in our data. Figure 
[9] plots our new contrast and mass limits as a function of 
projected stellocentric dis tance; we adopt the Hipparcos 
system distance of 157 pc ( van Leeuwen|2007 ). We reach 
a sensitivity of 20 Mj„p beyond 100 AU and < 5 Mj„p 
beyond 175 AU. In the east half of the system, our FOV 
extended to 800 AU, while on the west side, the FOV 
ended at ~ 200 AU. An equally luminous binary com- 
panion to "B" (i.e. "Bb") is ruled out at > 15 AU sepa- 
ration from "B." 

4. DISCUSSION 
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J, H , and K flux relative to the L' and M' flux. From 
this we conclude that hot dust must be present within 
several BD radii away from the BD's surface, consistent 
with what we expect for an actively accreting disk. 

Adding additional weight to the disk hypothesis is 
the 24 /im excess. As described in Section |3.1.4[ the 
GSC 06 214 system has an unresolved 1.6 mJy excess at 



24 ^.m. Luhman & Mamajek (2012 1 recently compiled 
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Figure 9. LMIRCam V Scr contrast curve for HIP 78530 (black 
line). Red diamond is HIP 78530B. Dashed lines are apparent mag- 
nitudes for objects of various masses according to 5 Myr DUSTY 
models. The dot-dashed line denotes the extent of the FOV on the 
west half of the system. 

4.1. GSC 06214B Disk Properties 

GSC 06214B was already known to exhibit an L' excess 
as we ll as Pa/3 emission, s uggesting a circum(sub) stellar 



disk ( [Bowler et al ||2011 | 

pothesis using the full co 



We wished to test this hy- 
lection of 1 — 5 fj,m and 24 fj,m 
photometry. The shape of the 1 — 5 /im excess SED is 
strongly dependent on the temperature (location) of the 
inner edge of the disk. An accreting dust disk would 
likely be truncated at the dust sublimation temperature 
or by magnetic interactions near the c orotation radius 



Ks — [24] and [4.5J — [24J colors for other low-mass USco 
objects. The ratios of the 24 fj,m excess to the K and M' 
fluxes of GSC 06214B fall within the range that the au- 
thors derived for other disk-bearing late M-type objects 
in the association. Although the primary is shown to 
have Ha emission, presumably due to an accretion disk, 
the data are consistent with both the near- and mid-IR 
broadband excesses originating from a single, continuous 
disk around the BD. 

4.2. HIP 78530B Spectral Type and Membership 

As discussed in Section [3A2l HIP 78530B's broad- 
band colors do not match the M8 spectral class indicated 
by the shape of its J/H/K spectra. The L2010 young 
photosphere color calibration suggests it is an M3 ± 2 
spectral type (Figures |4| andjsl , and the observed colors 
also more closely match M5V colors than M8V colors 
(SPEX standards Gl 388 and LP 412-31 respectively). 
Using the temperature scale calibrated by [Luhman et al.] 
(20031, M3 corresponds to an effective temperature of 
SMT- 3400 K. The BCAH 5 Myr evolutionary models 

0.4 Mo, 



(Baraffe et al. 



(Shu et al. 



1994 



2003), both of which 



Muzerolle et al. 
are a few times the BD radius. 

To estimate the location of the disk inner edge, we cre- 
ated simple model rings of blackbody grains at a single 
temperature. A disk with an inner edge at a temperature 
of 1550 K, the sublimation temperature of silicates, re- 
produced the shape of the observed J through M' excess. 
Decreasing the temperature to 1000 K under-produced 



1998) predict a mass of 0.18 
depending on whether the 3350 K effective temperature, 
the H — K color, or the H — L' color is used. Regard- 
less of the particular assumptions, these models place the 
object above the hydrogen burning limit (^ 75 Mj). 

Conversely, the apparent magnitude of this object is 
consistent with other M8 brown dwarfs in Upper Sco 
(Figure |8]) as well as the DUSTY and BCAH predictions 
for a 5 -T.0 Myr, 20 Mjup BD. The 13 pc uncertainty in 
distance to HIP 78530A cannot account for the faintness 
of the object compared to other USco M3 BDs. New pho- 
tometry at M' would help determine whether the object 



Environments of USco Brown Dwarf Companions 



11 



is blue at all wavelengths or merely faint at K and L' . 

Although the primary has a Hipparcos parallax which 
is consistent with the distance to USco, it is possible 
that "B" is a background ob ject. In the discovery pa- 
per, [Lafreniere et al. (|2011 ) found agreement between 
the strength oi gravity-sensitive NIR spectral features in 
HIP 78530B and those in another USco M8 BD. Addi- 
tionally, the authors measured common proper motion 
at the 6(T level over a 2 year baseline. Unfortunately 
our astrometric calibration is not sufficiently precise to 
confirm or deny common proper motion. An additional 
epoch of observation would help to cement this object's 
status as a bound companion. 

5. DISKS AS CONSTRAINTS ON SCATTERING EVENTS 

The presence of an accreting disk around a low-mass 
companion at a wide separation from its primary star has 
the potential to limit the possible formation histories of 
the companion. To see why this is the case, we consider 
the fate of a gas disk orbiting a brown dwarf of mass 
M, which is in turn orbiting a star of mass M*. An 
encounter or series of encounters with a perturber of mass 
Mp > M, also orbiting the primary star, can scatter the 
brown dwarf to wide separations, while typically leaving 
the more massive perturber on an orbit near the primary. 

Each encounter truncates the disk at a distance, 
r, from the brown dwarf. This truncation distance 
is roughly the location where the tidal force dur- 
ing the encounter generates a velocity kick Aw ~ 
{GMp/b^){r/b){2b/v^nc) equal to the escape velocity 
from the brown dwarf, Uosc — (2G'Af/r)^/^, where G is 
the gravitational constant, and b and Vonc are the dis- 
tance of closest approach and the relative velocity of the 
two bodies during the encounter, respectively. We con- 
sider a fiducial encounter for which b is of order the Hill 
radius of the perturber, b ~ Rh,p ~ a{Mp/3M^y/'^ and 
Venc = ^b, is the relative velocity of two circular Keple- 
rian orbits with semi-major axes that differ by b. Here, a 
is the distance from the host star at which the encounter 
occurs, and n = iGNh/a^f/'^ is the angular velocity 
of a circular orbit at that distance. These choices are 
intended to maximize the expected disk truncation dur- 
ing scattering, without invoking unusually close or slow 
encounters, which we expect to be rare. Note that ar- 
bitrarily extreme encounters are disallowed because the 
BD cannot be ejected from the system. Setting Aw = Vesc 
yields a truncation radius of r ^ 0.5Rh, where Rh is the 
brown dwarf's Hill radius. The disk truncation radius 
from scattering is thus predicted to be comparable to 
the radius generated merely by long-t erm tidal interac 
tions with the host star, ~ OARh (e.g. Martin & Lubow 



2011 1. To solidify interpretation of future observations. 



we suggest that numerical simulations of scattering histo- 
ries producing wide-separation BD companions should be 
mined to determine the typical properties of their most 
disruptive close encounters. 

Although scattering does not significantly strip disk 
material from the BD, the disk can still be used to con- 
strain the BD's scattering history. This is because the 
disk size, ~ OARh-, is directly proportional to a, the 
semi-major axis at which the BD formed. We next esti- 
mate the range of formation locations which would pro- 
duce disks large enough to survive for 5 Myr, the mini- 
mum age of USco. 



Observations of protoplanetary disks at separations of 
tens of AU yield surf ace densitv profiles o f order S = 
10^ g/cm^ (r/AU)~^ (Andrews et al. 20091, comparable 
to the minimum-mass solar nebula. L)isks around brown 
dwarfs are likely less massive, so this profile places an 
upper limit on the disk lifetime. If such profiles extend 
inward to smaller separations from the BD, a typical disk 
contains a mass of Md ~ 7 x lO~'*M0(r'/AU) within a 
radius r. The size of the inner clearing of a n accreting 
disk is negligible, as we found in Section [4.1| Assuming 
that a typical low-mass brown dwarf has an accretion 
rate comparable to G SC 06214B, Af = 2 x lO-^Me/yr 
(Bowler et al. 20111, it would accrete at least 10"'* 
of material over 5 Myr, the estimated mass interior to 
r — 0.15 AU. A higher accretion rate of lO~^°M0/yr 
would process 5 x lO~'*Af0 over 5 Myr, equal to the 
estimated mass within r = 0.7 AU. 

These estimates of the minimum disk size can be 
used to constrain the formation location of the BD. 
A truncation radius of 0.4a(Af/3M*)*/"^ corresponds to 
0.06a — 0.07a for our three BDs, using the masses listed 
in Tables [2j [3] and |4j We estimate that insufficient mass 
would remain in the current disk orbiting GSC 062 14B 
if the brown dwarf were scattered outward by a per- 
turber or perturbers from a stellocentric distance of a < 
2 — 10 AU. Current limits on additional companions in 
the GSC 06214 system can rule out 35Mj„p ~ 2M^ ob- 



jects beyond 10 AU stellocentric distance (Section 3.1.5 1, 
disfavoring, though note entirely ruling out, the scatter- 
ing hypothesis. Limits on the minimum semi-major axis 
at formation of HIP 78530B and IRXS 1609B, based on 
their lack of accretion signatures, are similar, presuming 
they had accretion rates similar to that of GSC 06214B. 
Additional objects more than twice as massive as the 
known BDs are ruled out beyon d 35 AU for IRXS 1609 
and beyon d 75 AU for HIP 78530 ( |Lafreniere et al.| ( |2010 1 
and 3.3.41. These weaker constraints cannot be used to 
rule out a scattering origin for either object. Given our 
estimates, tighter observational constraints on the pres- 
ence of possible perturbers, coupled with more detailed 
modeling of disk truncation during ejection by scattering 
will provide significant limits on the scattering histories 
of wide-separation, low-mass companions. 

6. SUMMARY AND FUTURE WORK 

We present a 3-5 fxm LBT/MMT imaging study of 
three low-mass (ratio) companions to Upper Scorpius 
stars: GSC 06214B, IRXS 1609B, and HIP 78530B. 
These three systems constitute some of the lowest mass- 
ratio {q ~ 1%) non-planetary systems discovered to date. 
However, though the companions have similar separa- 
tions and mass ratios, they probe two very different types 
of binary systems: GSC' 06214 and IRXS 1609 are K- 
type/L-type pairs, while HIP 78530 is a B-type/M-type 
pair. Combined with MIPS-24 /zm photometry, we use 
our new high-resolution near/thermal IR data to iden- 
tify evidence for broadband disk emission from small 
dust grains around both the companions and the primary 
stars. Our study yields the following major results: 

• We confirm the presence of a circum(sub)stellar 
disk around the LO-type BD GSC 06214B. Fur- 
thermore, we find the best fit temperature of the 
disk inner edge is at the dust sublimation tempera- 
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ture, consistent with active accretion. A moderate 
24 iim excess in the system is plausibly explained 
by the contribution from the BD's disk, although 
some material must still be present around the pri- 
mary in order to explain its Ha emission. The ex- 
istence of a disk around GSC 06214B suggests, but 
does not conclusively prove, it has not undergone 
a scattering event. 

• We find evidence of a weak 24 /im excess in the 
IRXS 1609 system, but the < 4 fim photome- 
try of the primary and low-mass companion are 
consistent with purely photospheric emission. The 
24 fj,m excess is unresolved, so we cannot determine 
whether "A," "B," or both host a warm disk. 

• In contrast, we find no evidence of 1 — 4 yum or 
24 fj,m excesses in the HIP 78530 system. How- 
ever, we find a discrepancy between the spectral 
types indicated by the 1—4 fim broadband colors 
and by the J/H/K spectra. Although the spec- 
tra and apparent magnitude are consistent with an 
USco M8-type BD, the broadband colors instead 
match a young M3-type low-mass star. Future 
work corroborating common proper motion in the 
system is necessary to decisively rule out the possi- 
bility that "B" is a background object. We present 
improved constraints on the presence of additional 
low-mass companions, ruling out equally luminous 
companions at projected separations from the pri- 
mary > 100 AU and reaching a background limit 
of < 5 Mjup beyond 175 AU . 

Although our sample size is too small for robust statis- 
tics, the detection of a massive, hot disk in one of three 



tion determined for single BDs in USco 1 


Riaz et al.|2012 


Scholz et al.||2007 


Luhman & Mamajek 


201:^). ^'urther- 



mation mechanisms for low-mass, widely separated com- 
panions remains unclear, this work does provide some 

constraints on the formation of our USco targets^ 

In particu lar, we corroborate the findings of |Bowler| 
eFaLlpOTT] ), who identified a disk around GSC 06214B 
from accretion-driven Pa/3 emission. Our study fur- 
ther confirms that this disk is dusty, producing a strong 
broadband IR excess from grains near the dust subli- 
mation radius, much like for many disks around stars 



al. 2006 2009 Luhman & Ma- 



in USco ([Carpenter et 

majek 20l2p . I'he presence ot a dusty, accreting disk 
around the BD, combined with the constraint that no 
object more than twice as massive as "B" exists beyond 
10 AU stellocentric distance, casts doubt on whether this 
companion could have be en formed via a p lanet-planet 
scattering event (i.e. as in [Veras et al.|[2009l ). 

We do not detect a near-lR broadband excess around 
IRXS 1609B. (The fact that the L4-type IRXS 1609B 
does not have a hot disk, while the more massive LO- 
type GSC 06214B does, need not be contradictory; it is 
within the scatter in the inner disk dissipation timescale.) 
We do, however, find the system as a whole exhibits a 
24 fim excess, presumably from warm dust. If some of 
this emission originates from a large disk surrounding the 
BD companion, then a scattering origin for this BD may 
also be disfavored. 



Conversely, a lack of 24 excess around the more 
massive HIP 78530B does not necessarily imply a past 
scattering event, as depletion of warm dust from passive 
disk evolution alone is more rapid for intermediate- mass 
than for low-mass BDs ( |Luhman fc Mamajekf2012l ). 

Resolved observations with ALMA and/ or IN (J MIC 
(the 8— 13 fim channel of LBTI) would be able to further 
constrain the disk properties of both the GSC 06214 pri- 
mary and BD and could clarify whether the IRXS 1609 
system's warm 24 //to excess emission originates solely 
from material around the primary or around the low- 
mass companion as well. ALMA observations could in 
principle identify even colder dust emission from outer 
disk regions and provide a direct measurement of disk 
mass. They might also be used to differentiate between 
fully depleted disks and those with central clearings. A 
comparison of the incidence, masses, and morphologies 
derived in this way for disks around a larger sample of 
single and companion BDs might discriminate between 
in situ formation and scattering events. Finally, com- 
paring cluster mass functions with age may identify a 
time-dependent low-mass BD population consistent with 
some fraction of low- mass BDs having a scattering origin. 
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